The Astrophysical Journal, 2005 October, Iastro- ph/0506608 

Preprint typeset using PJTpX style emulateapj v. 7/8/03 



THE X-RAY SPECTRAL PROPERTIES OF SCUBA GALAXIES 
D.M. Alexander, 1 F.E. Bauer, 1 S.C. Chapman, 2 I. Smail, 3 A.W. Blain, 2 W.N. Brandt, 4 and R.J. Ivison 5 ' 6 

Institute of Astronomy, Madingley Road, Cambridge CB3 OHA, UK 
California Institute of Technology, Pasadena, CA 91125, USA 
3 Institute for Computational Cosmology, University of Durham, South Road, Durham DH1 3LE, UK 
4 Department of Astronomy and Astrophysics, Pennsylvania State University, 525 Davey Laboratory, University Park, PA 16802, USA 
5 Astronomy Technology Centre, Royal Observatory, Blackford Hill, Edinburgh EH9 3HJ, UK and 
6 Institute for Astronomy, University of Edinburgh, Blackford Hill, Edinburgh EH9 3HJ, UK 
Received 2005 Feb 22; accepted 2005 June 23 

ABSTRACT 

Deep SCUBA surveys have uncovered a large population of massive submillimeter emitting galaxies (SMGs; 
/850^m ~ 4 mJy) at z ^ 1 . Although it is generally believed that these galaxies host intense star-formation 
activity, there is growing evidence that a substantial fraction also harbor an Active Galactic Nucleus [AGN; 
i.e., an accreting super-massive black hole (SMBH)]. We present here possibly the strongest evidence for this 
viewpoint to date: the combination of ultra-deep X-ray observations (the 2 Ms Chandra Deep Field-North) and 
deep Keck spectroscopic data of SMGs with radio counterparts. We find that the majority (ss 75%) of these 
radio-selected spectroscopically identified SMGs host AGN activity; the other « 25% have X-ray properties 
consistent with star formation (X-ray derived star-formation rates of « 1300-2700 M yr~'). The AGNs have 
properties generally consistent with those of nearby luminous AGNs (Ffa 1 .8 ± 0.5, A^h ~ 10 20 — 10 24 cm~ 2 , and 
Lx ~ 10 43 — 10 44 ' 5 erg s _1 ) and the majority (w 80%) are heavily obscured (Nh ^ 10 23 cm -2 ). We construct 
composite rest-frame 2-20 keV spectra for three different obscuration classes (A^h < 10 23 cm~ 2 , A^h = 1— 
5 x 10 23 cm~ 2 , and A^h > 5 x 10 cm~ 2 ) which reveal features not seen in the individual X-ray spectra. An 
k 1 keV equivalent width Fe Koc emission line is seen in the composite X-ray spectrum of the most heavily 
obscured AGNs, suggesting Compton-thick or near Compton-thick absorption. Even taking into account the 
effects of absorption, we find that the average X-ray to far-infrared luminosity ratio of the AGN-classified 
SMGs (j^= 0.004) is approximately one order of magnitude below that found for typical quasars. This result 

suggests that intense star-formation activity (of order « 1000 M Q yr~') dominates the bolometric output of 
these SMGs. However, we also explore the possibility that the X-ray to far-infrared luminosity ratio of the 
AGN components is intrinsically less than that found for typical quasars and postulate that some SMGs may be 
AGN dominated. We investigate the implications of our results for the growth of massive black holes, discuss 
the prospects for deeper X-ray observations, and explore the scientific potential offered by the next generation 
of X-ray observatories. 

Subject headings: infrared: galaxies — X-rays: galaxies — galaxies: active — galaxies: starburst 



1. INTRODUCTION 

Deep SCUBA surveys have uncovered a large population 
of submillimeter (submm; A, = 300-1000 /jm) emitting galax- 
ies (SMGs; <, 1000 sources deg~ 2 at /sso^m ^ 4 mJy; e.g., 
Smail, Ivison, & Blain 1997; Barger, Cowie & Sanders 1999; 
Blain etal. 1999, 2002; Eales etal. 1999; Cowie etal.2002; 
Scott etal.2002; Borys etal. 2003; Webb etal. 2003a). The 
majority of these SMGs are faint at all other wavelengths, hin- 
dering source identification studies (e.g., Hughes etal. 1998; 
Ivison et al. 2000; Smail et al. 2002; Webb et al. 2003b). How- 
ever, due to a considerable amount of intensive multi- 
wavelength follow-up effort, it is becoming clear that al- 
most all are dust-enshrouded luminous galaxies at z > 1 
(e.g., Ivison etal. 1998, 2000; Barger etal. 1999; Chapman 
etal. 2003a, 2005; Simpson etal. 2004; Pope etal. 2005). 
With estimated bolometric luminosities of order 10 13 L 
(Blain etal. 2004a; Chapman etal. 2005), SMGs outnum- 
ber comparably luminous galaxies in the local Universe by 
several orders of magnitude (e.g., Cowie etal. 2004). Dy- 
namical mass and galaxy clustering studies of spectroscop- 
ically identified SMGs have indicated that they are likely 
to be massive galaxies (e.g., Frayer etal. 1998, 1999; Gen- 
zel etal. 2003; Neri etal. 2003; Blain etal. 2004b; Swin- 



bank etal. 2004; Greve etal. 2005); these SMGs are w 5 
times more massive than the coeval optically selected field 
galaxies (e.g., Erb etal. 2003; Steidel etal. 2004; Swinbank 
etal. 2004). Making reasonable assumptions about the dura- 
tion of the submm-bright phase in these systems, these re- 
sults provide evidence that SMGs are the progenitors of lo- 
cal ^M* early-type galaxies (e.g., Lilly etal. 1999; Scott 
et al. 2002; Smail et al. 2004; Chapman et al. 2005). 

Central to the study of SMGs is the physical origin of their 
extreme luminosities [i.e., starburst or Active Galactic Nu- 
cleus (AGN) activity]. If these sources are shown to be ultra- 
luminous starburst galaxies then their derived star-formation 
rates suggest a substantial increase in star-formation activity 
at z > 1 (e.g., Blain et al. 1999; Smail et al. 2002). Conversely, 
if these sources are shown to be powerful AGNs then they out- 
number comparably luminous optical quasars by « 1-2 orders 
of magnitude. The apparent association of a few SMGs with 
quasars and their comparable comoving space densities (when 
corrected for probable source lifetimes) suggest an evolu- 
tionary connection between these two populations (e.g., Page 
etal. 2001, 2004; Croom et al. 2004; Stevens et al. 2003, 2004; 
Chapman et al. 2005). Given this evidence it is plausible that 
SMGs represent a rapid black-hole and stellar growth stage 
prior to a luminous quasar phase (e.g., Archibald etal. 2002; 
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Almaini 2003; Page et al. 2004; Alexander et al. 2005a). 

Arguably the most direct indication of AGN activity is the 
detection of luminous hard X-ray emission (i.e., > 2 keV). 
Hard X-ray emission appears to be a universal property of 
AGNs, giving a direct "window" on the emission regions 
closest to the SMBH (e.g., Mushotzky, Done, & Pounds 
1993), and it can provide a secure AGN identification in 
sources where the optical signatures and counterparts are 
weak or even non existent (e.g., Alexander etal. 2001; Co- 
mastri etal. 2002). Hard X-ray emission is relatively insensi- 
tive to obscuration (at least for sources that are Compton thin; 
i.e., A^h 1 .5 x 10 24 cm~ 2 ) and any hard X-ray emission from 
star formation in the host galaxy is often insignificant when 
compared to that produced by the AGN. 

The first cross-correlation studies of moderately deep 
(« 20-200 ks) Chandra surveys with SCUBA surveys 
yielded little overlap between the X-ray and submm detected 
source populations (^10-20%; e.g., Fabian etal. 2000; 
Bautz etal. 2000; Hornschemeier etal. 2000; Severgnini 
etal. 2000). 1 To first order these results suggested that bolo- 
metrically dominant AGNs can only be present in typical 
SMGs if they are Compton thick. Later studies with the Chan- 
dra Deep Field-North (CDF-N; Brandt etal. 2001, Alexan- 
der et al. 2003a) survey showed that a significant fraction (up- 
wards of « 30-50% when the Borys et al. 2003 SCUBA map 
is used) of SMGs are X-ray detected (Barger etal. 2001b; 
Alexander et al. 2003b; Borys et al. 2004). Direct X-ray spec- 
tral analyses of the five AGN-classified SMGs in Alexander 
etal. (2003b) indicated that the AGNs were heavily obscured 
but only moderately luminous at X-ray energies (Lx ~ 0.3- 
1x10 erg s _1 when corrected for the effects of absorp- 
tion). A comparison of the X-ray-to-submm spectral slopes of 
these SMGs to that of NGC 6240 (a nearby luminous galaxy 
with an obscured AGN) suggested that the AGNs typically 
contributed only a few percent of the bolometric luminosity. 
However, the small sample size and lack of spectroscopic red- 
shifts (only one source had a spectroscopic redshift) prevented 
more quantitative conclusions. 

The deep optical spectroscopic work of Chapman 
etal. (2003a, 2005) has recently provided spectroscopic red- 
shifts for 73 radio-identified SMGs, a significant increase 
in sample size over previous studies (e.g., Ivison etal. 1998, 
2000; Barger et al. 1999; Ledlow et al. 2002; Smail et al. 2003; 
Simpson etal. 2004). The 2 Ms CDF-N field was one of 
the regions targeted for this intensive spectroscopic follow 
up. The combination of deep optical spectroscopic data and 
ultra-deep X-ray observations provides powerful constraints 
on AGNs in SMGs. In particular, spectroscopic redshifts im- 
prove the accuracy of the X-ray spectral analyses over those of 
Alexander et al. (2003b) through the identification of discrete 
X-ray spectral features (e.g., Fe Ka emission) and the deter- 
mination of any intrinsic absorption, which is a strong func- 
tion of redshift in a given X-ray band [Nh,z ~ Nn, z o(^ +z) 2 ' 6 ]- 
These improvements promise the most accurate determina- 
tion of the AGN contribution to the bolometric output of these 
SMGs to date. In this paper we investigate the X-ray proper- 
ties of the radio-selected spectroscopically identified SMGs 
in the CDF-N field and predict the AGN contribution to the 
bolometric luminosity of these SMGs. The Galactic column 

1 Further cross-correlation studies with moderately deep X-ray observa- 
tions have revealed some overlap (e.g., Barger etal. 2001a; Ivison etal. 2002; 
Almaini et al. 2003; Waskett et al. 2003). This was probably due to the larger 
areal coverage of the SCUBA observations in these studies. 



density toward the CDF-N field is (1.3 ±0.4) x 10 20 cm" 2 
(Lockman 2004), and H = 65 km s -1 Mpc~\ ii M = 3, and 
i2 A = I are adopted throughout. 



2. OBSERVATIONS AND BASIC ANALYSIS 

2.1. SCUBA Galaxy Sample 

Our SCUBA galaxy sample includes 20 SMGs that have 
been spectroscopically identified by Chapman et al. (2005) 
and lie in the region of the 2 Ms CDF-N observation (Alexan- 
der et al. 2003a). The optical counterparts for these sources 
were chosen based on the positional association with a radio 
source [identified by either VLA (Richards 2000) or MER- 
LIN (Chapman et al. 2004a) observations]. The basic proper- 
ties of the SMGs are given in Table 1 . The source redshifts 
(z = 0.555-2.914 with a median of z = 2.0±0.6) are typical 
of the spectroscopically identified SMG population (median 
redshift of z « 2.2; Chapman etal. 2005). The w 35-50% 
of sources without radio counterparts (e.g., Ivison et al. 2002; 
Chapman et al. 2003b; Borys et al. 2004) could not be targeted 
with optical spectroscopy. However, the range of 850/jm 
flux densities for the submm sample (/850/jm ~ 4-12 mJy, 
with a mean of 6.6 ± 2.2 mJy) is typical of SMGs detected 
in blank-field SCUBA surveys (e.g., Scott etal. 2002; Webb 
etal. 2003a). 

The potential selection effects introduced by the radio and 
spectroscopic-identification of this sample, compared to a 
purely submm-flux-limited sample of SMGs, are discussed at 
length in Chapman etal. (2005) and Blain etal. (2004). The 
need for both a faint radio counterpart and identifiable spec- 
tral features to aid in redshift measurement may result in this 
sample having a higher incidence of AGN activity than the 
SMG population as a whole; however, the incompleteness in 
our surveys due to these two selection criteria can be equally 
well explained by the likely temperature and redshift distribu- 
tion of the SMG population (Chapman et al. 2005). Neverthe- 
less, we caution the reader that the sample analysed here may 
not be completely representative of the entire /sso^m ^ 4 mJy 
SMG population. 

Many of the SMGs in our sample were radio-selected 
sources specifically targeted with SCUBA observations; see 
the table in Alexander etal. (2005a) for the observation 
modes. While this could potentially cause an additional AGN 
bias over a purely radio-detected SMG sample, we do not 
find evidence for a strong bias in our sample (see Alexander 
etal. 2005a for the results of a two-sided Fisher's exact test). 

2.2. Chandra Counterparts to the SCUBA Galaxy Sample 

The 2 Ms CDF-N observations were centered on the opti- 
cal Hubble Deep Field-North (HDF-N; Williams etal. 1996) 
region and cover w 448 arcmin 2 (Alexander et al. 2003a). 
These observations provide the deepest view of the Uni- 
verse in the 0.5-8.0 keV band; the aim-point sensitivities 
are « 7.1 x 10~ 17 erg cm~ 2 s _1 at 0.5-8.0 keV (full band), 
« 2.5 x 10~ 17 erg cm~ 2 s _1 at 0.5-2.0 keV (soft band), and 
« 1.4 x 10~ 16 erg cm" 2 s" 1 at 2-8 keV (hard band). At the 
median redshift of the SMGs, the observed 0.5-8.0 keV band 
corresponds to rest-frame energies of 1 .5-24 keV. Such high 
X-ray energies allow very high column densities to be pene- 
trated (e.g., at 20 keV, < 55% of the direct X-ray emission 
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is absorbed with column densities of A^h 10 24 cm~ 2 ; see 
Appendix B in Deluit & Courvoisier 2003). 

Using a 1.5" search radius, 16 of the 20 SMGs were found 
to have Chandra counterparts in the main Chandra cata- 
log of Alexander etal. (2003a). Given the low surface den- 
sity of SMGs, we also searched for X-ray counterparts us- 
ing the complete supplementary Chandra catalog of Alexan- 
der et al. (2003a) and found one further match; the probability 
of this match being spurious is < 1% (see §2.4 of Alexander 
et al. 2003b). In total 17 (85±^%) of the 20 SMGs in our 
sample have a Chandra counterpart (see Table l). 2 We have 
calculated rest-frame 0.5-8.0 keV luminosities, and 1.4 GHz 
luminosity densities following equations 1 and 2 in Alexander 
et al. (2003b), assuming an X-ray spectral slope of F = 1 .8 and 
a radio spectral slope of a = 0.8, respectively. Far-infrared 
(far-IR; X = 40-120 /jm) luminosities were calculated from the 
rest-frame 1 .4 GHz luminosity densities under the assumption 
of the radio-to-far-IR correlation with q = 2.35 (e.g., Helou, 
Soifer, & Rowan-Robinson 1985); the presence of an AGN 
component to the radio emission will lead to an overpredic- 
tion of the far-IR luminosity. See Table 1 . 

The three X-ray undetected SMGs lie in sensitive regions 
of the CDF-N field and have 3 O flux limits below the 
fluxes of the X-ray detected SMGs; see Table 1. Their red- 
shifts (z = 2.0-2.1) are consistent with the median redshift 
of the SMG sample. Stacking the X-ray data of the individ- 
ually undetected SMGs following the procedure of Lehmer 
et al. (2005) yields marginally significant detections in the soft 
and full bands (2.5a, soft band; 2.7c, full band; S. Immler, 
private communication), suggesting that these sources lie just 
below the individual source detection limit. The correspond- 
ing average X-ray constraints are ss 9.0 x 10~ 17 erg cm~ 2 s _1 
(full band), w 1.7 x 10~ 17 erg cm~ 2 s _1 (soft band) and 
< 1.6 x 10~ 16 erg cm~ 2 s _1 (hard band). The full-band con- 
straint corresponds to Lq 5_8 o ~2 x 10 42 erg s~ ! at z = 2 
for r = 1 .8 and is consistent with that expected from luminous 
star formation (see §3.2 and §4.2). 

2.3. X-ray Spectra: Extraction and Fitting 

X-ray spectra were generated for all of the SMGs with full- 
band fluxes > 3 x 10~ 16 erg cm~ 2 s -1 , which typically cor- 
responds to > 50 X-ray counts. To account for the range of 
roll angles and aim points in the 20 separate observations that 
comprise the 2 Ms CDF-N, the X-ray spectra were generated 
using the ACIS source extraction code (ACIS Extract) de- 
scribed in Broos et al. (2002). 3 Briefly, for each source this 
code extracts the counts from each of the observations tak- 
ing into account the changing shape and size of the PSF with 
off-axis angle as given in the Chandra X-ray Center (CXC) 
PSF library. 4 A local background is extracted after all of the 
sources are excluded from the X-ray event file, and the spec- 
tra and response matrices are summed using standard FTOOLS 
routines (Blackburn 1995). See F. E. Bauer et al. (in prepara- 
tion) for further information and the X-ray spectral properties 
of the sources detected in the Chandra deep fields. 

Given the limited counting statistics of the X-ray sources, 
we performed all of the X-ray spectral analyses using the C- 
statistic (Cash 1979). One advantage of using the C-statistic 

2 All errors are taken from Tables 1 and 2 of Gehrels (1986) and correspond 
to the 1<J level; these were calculated assuming Poisson statistics. 

3 ACIS EXTRACT is a part of the TARA software package and can be ac- 
cessed from http://www.astro.psu.edu/xray/docs/TARA/ae_users_guide.html 

4 See http://asc.harvard.edU/ciao2.2/documents_dictionary.html#psf 
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FIG. 1 . — Submm flux density versus full-band flux for the SMGs (filled 
circles) and optically classified quasars with X-ray and submm constraints 
(open triangles). The quasar data are taken from Page etal. (2001), Vignali 
etal. (2001), and Isaak etal. (2002). The dashed line indicates the expected 
X-ray fluxes and submm flux densities for a quasar with the same properties 
as 3C273; the submm-to-X-ray spectral slope is independent of redshift (see 
Figure 2 of Fabian etal. 2000). Our spectroscopically identified SMGs are 
up to two orders of magnitude fainter in the X-ray band than the optically 
classified quasars for a given submm flux density. 



is that the data can be fitted without binning, making it ide- 
ally suited to low-count sources (e.g., Nousek & Shue 1989). 
A disadvantage of using the C-statistic is that it is not pos- 
sible to determine rigorously if one model provides a sta- 
tistically preferable fit than another; however, it is possible 
to inspect the fit residuals and perform Monte-Carlo analy- 
ses to distinguish between different models. Hence in our 
analyses we have focused on empirically motivated models 
which have proven to provide a robust characterisation of 
the properties of well-studied AGNs. We used the latest 
version of XSPEC (vll.3.1; Arnaud 1996, 2002) for all of 
the model fitting, which allows the C-statistic to be used on 
background-subtracted data; we carried out several checks of 
the background-subtraction method to verify that no spurious 
residual features were present in the background-subtracted 
data. We fitted the X-ray spectra of the SMGs using a power 
law model (with Galactic absorption) in the rest-frame 2- 
10 keV and 5-20 keV bands; see Table 2. The constraints 
on the five sources with ^100 counts are poor. However, 
our analyses have been designed to maximize the useful con- 
straints for even the weakest X-ray sources. The fit parameter 
uncertainties are quoted at the 90% confidence level for one 
parameter of interest (Avni 1976). We provide details of our 
analysis techniques below but defer the discussion of our re- 
sults to §3. 

In order to interpret the X-ray spectral fitting results we 
have developed an empirically motivated AGN model and 
constructed a simple rest-frame power-law diagnostic dia- 
gram. The advantage of this approach is that we only need 
to perform simple rest-frame power law fits to the sources, al- 
lowing us to determine the characteristics of sources with few 



4 



ALEXANDER ET AL. 



10 



(a) 
AGN 

SF 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 
Classification #1 

Obscured Unobscured 
HMXBs 
ULXs/LMXBs 




-1.5 -1 -0.5 0.5 1 1.5 2 
Effective Photon Index (r) 



2.5 




10 23 10 24 10 25 

Rest-frame 1.4GHz luminosity (W Hz- 1 ) 



FIG. 2. — Classification scheme for the X-ray emission from the SMGs. (a) Effective photon-index histogram. The approximate observed X-ray spectral slopes 
for a variety of different source types [split broadly into AGN and star formation (SF)] are shown (e.g., Nandra & Pounds 1994; Maiolino etal. 1998; Colbert 
etal. 2004). The only sources that produce extremely fiat or inverted X-ray spectral slopes (r < 0.5) are obscured AGNs. The typical uncertainties in the X-ray 
spectral slopes for the plotted sources are AT" w 0.3; see Table 1. (b) Rest-frame 0.5-8.0 keV luminosity versus 1.4 GHz luminosity density. The sources classified 
as obscured AGNs from Figure 2a are indicated with crosses; the X-ray luminosities have not been corrected for the effect of absorption. The dotted line shows 
the X-ray-radio relationship for star-forming galaxies whose X-ray emission is dominated by HMXBs (Persic et al. 2004); this relationship is converted to the 
0.5-8.0 keV band from the 2-10 keV band assuming T = 1.8. The X-ray emission from 15 (w 75%) of the 20 SMGs is dominated by AGN activity. We caution 
the reader that due to selection biases this does not directly indicate a 75% AGN fraction in the bright SMG population; see §3.2. 



X-ray counts. Our adopted AGN model includes a power-law 
component (which can be absorbed by neutral material; A^h), a 
neutral reflection component (the PEXRAV model in XSPEC), 
a scattered component (Sc) of ionised gas, and a Fe Ka emis- 
sion line at 6.4 keV. The XPSEC model components for our 
adopted model are (POW+PEXRAV+ZGA)+POW*ZWABS. 

We used the XSPEC command fakeJt to produce the model 
tracks for a range of X-ray spectral slopes (r = 1.3-2.3) and 
absorbing column densities (A^h = 10 -2x 10 24 cm -2 ). The 
ratio of the direct to reflected emission (i.e., reLrefl in the 
PEXRAV model) was fixed to unity, in agreement with the av- 
erage value found for local AGNs (e.g., Risaliti etal. 2002; 
Deluit & Courvoisier 2003; Malizia etal. 2003). We set the 
other parameters of the neutral reflection component to be 
the same as those in Model 2c of Malizia etal. (2003). The 
ionised gas scattering fraction was set to be 1% of the in- 
trinsic power-law emission, and the equivalent width of the 
Fe Ka emission line was varied depending on the absorbing 
column density (W\ = 0.145 keV for Wh< 3 x 10 23 crrr 2 , 
W x = 0.37 keV for N H = 3-7 x 10 23 cm -2 , and W x = 1 .0 keV 
for Nr > 7 x 10 23 crrr 2 ; see Risaliti 2002 and Malizia 
etal. 2003). We also investigated higher column densities 
with a varying amount of ionised gas scattering fraction [Na ~ 
10 25 cm~ 2 (i.e., pure reflection) and Sc = 0-10%]. To take ac- 
count of the Chandra ACIS-I instrumental response we used 
the response matrix file (RMF) and ancillary response file 
(ARF) from CXOHDFN J 1 23629. 1 +62 1 045 . To allow for the 
construction of the model tracks up to rest-frame energies of 
20 keV we assumed a redshift of z — 2.0, consistent with the 
average redshift of the SMG sample. 



We also attempted to fit directly the individual X-ray spec- 
tra of the SMGs with the empirically motivated model out- 
lined above. However, we could not statistically constrain 
these models sufficiently to gain additional physical insight 
over that found using the power-law diagnostic diagram anal- 
yses. This was most likely due to a degeneracy in the number 
of model components and the generally poor photon statistics 
of the individual X-ray spectra. 

3. X-RAY PROPERTIES OF THE SCUBA GALAXY SAMPLE 

3.1. Basic Source Properties 

We show the submm flux density versus full-band flux of all 
of the SMGs in Figure 1 . As a comparison we also show the 
submm and X-ray properties of samples of moderate-to-high 
redshift quasars (z ps 1-5; Page et al. 2001 ; Vignali et al. 2001 ; 
Isaak et al. 2002). The average redshift and submm flux 
density of the Page etal. (2001) quasars (z = 2.1 ±0.4 and 
/850^m = 7.5 ±1.6 mJy) are consistent with those of our 
SMGs; however, the median X-ray flux of the SMGs is rj 30 
times fainter. The submm properties of quasars are proving 
to be important from the point of view of the processes of 
star formation at high redshift and providing clues to the rela- 
tionship between the formation of SMGs, quasars, and mas- 
sive galaxies (e.g., Stevens etal. 2003, 2004; Page et al. 2004). 
However, the X-ray properties of spectroscopically identified 
SMGs can only be investigated with ultra-deep X-ray obser- 
vations. 

3.2. X-ray Source Classification 

Due to the exceptional sensitivity of the 2 Ms CDF-N obser- 
vations it is possible to detect X-ray emission from star forma- 
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tion out to high redshift. For example, the soft-band flux limit 
at the aim point corresponds to Z4.5-6 kev > 8 x 10 41 erg s _1 
for a galaxy at z = 2. This luminosity is comparable to 
those for the most X-ray luminous starburst galaxies in the 
local Universe (e.g., NGC 3256; Moran etal. 1999; Lira 
et al. 2002), sources that are bolometrically less energetic than 
SMGs by w 1 order of magnitude. Therefore, the detection 
of X-ray emission from an SMG does not a priori indicate 
the presence of AGN activity. In order to classify the SMGs 
in our sample we have used a two-tiered classification pro- 
cedure. This classification procedure is conservative and may 
not identify all of the SMGs that host AGN activity. However, 
we can be confident that any AGN classifications are secure. 

The first method of source classification is based on the 
effective X-ray spectral slope; see Figure 2a. Star-forming 
galaxies and unobscured AGNs generally have comparatively 
steep X-ray spectral slopes (i.e., r ks 2; e.g., Kim, Fab- 
biano, & Trinchieri 1992a,b; Nandra & Pounds 1994; Ptak 
et al. 1999; George et al. 2000; Colbert et al. 2004). 5 Obscured 
AGNs are usually distinguishable from unobscured AGNs and 
star-forming galaxies by the presence of a flat X-ray spec- 
tral slope (T < 1) due to the energy-dependent photo-electric 
absorption of the X-ray emission (e.g., Maiolino etal. 1998; 
Risaliti, Maiolino, & Salvati 1999). However, since SMGs 
are potentially massive galaxies undergoing intense star for- 
mation, they could have a large population of high-mass X-ray 
binaries (HMXBs), which may have comparatively flat X-ray 
spectral slopes (e.g., neutron star HMXBs can have T « 0.5- 
1.0; Colbert etal. 2004). Eleven of the 17 X-ray detected 
SMGs have F < 0.5, unambiguously indicating the presence 
of a heavily obscured AGN. 

The second method of source classification focuses on the 
expected X-ray emission from star formation; see Figure 2b. 
Many studies have shown that the radio luminosity can be 
used to predict the X-ray luminosity in star-forming galaxies 
(e.g., Shapley, Fabbiano, & Eskridge 2001; Bauer etal. 2002; 
Ranalli etal. 2003; Grimm, Gilfanov, & Sunyaev 2003; Gil- 
fanov, Grimm, & Sunyaev 2004; Persic etal. 2004). The 
main assumption in these predictions is that the radio emis- 
sion is dominated by star formation; however, the presence 
of an AGN component to the radio emission will only lead 
to an overprediction of the contribution from star formation 
at X-ray energies. The normalisation and slope of these 
X-ray-radio relationships differ depending upon the relative 
contributions from low-mass and high-mass X-ray binaries 
(LMXBs; HMXBs; Gilfanov etal. 2004; Persic etal. 2004; 
see Figure 2b). Since the properties of SMGs indicate that 
they are undergoing intense star-formation activity, we only 
need to consider the contribution from HMXBs (see §4.2). 
The predicted X-ray emission from star formation in SMGs 
is large (up to Lx ~ 10 42 — 10 43 erg s _1 ). Even so, all of the 
sources classified as AGNs show an X-ray excess due to AGN 
activity; these X-ray excesses will be even greater once the 
X-ray luminosities are corrected for the effect of absorption. 
With this second method a further four sources show a clear 
X-ray excess over that predicted from star formation, indicat- 
ing that their X-ray emission is dominated by AGN activity. 
Four of the 10 AGN-classified SMGs with > 100 counts also 
show some evidence for X-ray variability (M. Paolillo, pri- 

5 Detailed X-ray studies of local star-forming galaxies have shown that 
they have more complex X-ray spectra than simple power-law emission (e.g., 
power-law emission and a very soft « 0.7 keV thermal component; Ptak 
etal. 1999); however, at the probable redshifts of our sources the very soft 
thermal component will be redshifted out of the Chandra energy band. 
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FIG. 3. — X-ray band ratio versus spectroscopic redshift for the X-ray clas- 
sified AGNs. The light and dark shaded regions show the range of expected 
band ratios for an unabsorbed and absorbed AGN, respectively. These regions 
were calculated assuming a T = 1.8 ±0.5 power law with differing amounts 
of absorption (as shown); these simple AGN models have been calculated 
using PIMMS Version 3. 2d. The error bars correspond to the uncertainties in 
the band ratio. This simple figure has limited diagnostic utility; however, it 
suggests that almost all of the AGNs are obscured. 



vate communication; see Table 1), a further signature of AGN 
activity; these constraints are consistent with those found for 
similarly bright AGNs in the Chandra Deep Field-South (Pao- 
lillo et al. 2004). None of the sources that have X-ray emission 
consistent with star formation show evidence for X-ray vari- 
ability, although the photon statistics are poor. 

Overall, the X-ray emission from 15 (75+ 2 g%) of the 20 
SMGs in the sample is AGN dominated. The other five 
SMGs are the faintest X-ray sources in the sample, and three 
sources are undetected in the X-ray band (see §2.2). The X- 
ray emission from these five sources is likely to be dominated 
by star-formation activity (see §4.2 for further constraints); 
however, heavily obscured or X-ray weak AGNs cannot be 
ruled out (e.g., compare with the similar z = 2.285 galaxy 
FSC 10214+4724; Alexander etal. 2005b). We caution the 
reader that, due to our sample selection and completeness, 
these results do not directly indicate a w 75% AGN fraction 
in the bright SMG population (the AGN fraction is probably 
> 38^Jq%; Alexander etal. 2005a, consistent with earlier esti- 
mates from the X-ray data; e.g., Alexander et al. 2003b; Borys 
et al. 2004; Wang, Cowie, & Barger 2004). 

3.3. X-ray Spectral Analyses 

The flat X-ray spectral slopes for the majority of the AGN- 
classified SMGs indicates the presence of absorption. Basic 
X-ray absorption constraints can be set using the X-ray band 
ratio (defined as the ratio of the hard-band to soft-band count 
rate); see Figure 3. This approach has limited diagnostic util- 
ity because AGNs often have more complex spectra than that 
of power-law emission with differing amounts of absorption. 
However, it suggests that the majority of the sources are heav- 
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FIG. 4. — Rest-frame 5-20 keV spectral slope versus rest-frame 2-10 keV 
spectral slope showing the X-ray properties of individual SMGs. The shaded 
region indicates where typical AGNs with T = 1.8 ±0.5 and Nh = 10 20 - 
2 x 10 24 cm~ 2 are expected to lie (individual values of Nn and T are indi- 
cated); see §2.3. The dot-dashed line shows where completely Compton- 
thick AGNs (i.e., JVh ~ 10 25 cm -2 ; a pure reflection spectrum with ionised 
gas scattering) with different amounts of scattering (as indicated) are ex- 
pected to lie; the assumed X-ray spectral slope is T = 1.8. This figure shows 
that the AGN-classified SMGs are well represented by the model, implying 
that the variations in their X-ray spectral properties are due to absorption. 



ily obscured (i.e., A^h ^ 10 23 cm -2 ). The intrinsic luminosity 
of the obscured AGNs could be considerably greater than the 
observed luminosity if the obscuration is high. 

In order to estimate the amount of absorption in each of 
the AGN-classified SMGs more accurately we have con- 
structed a rest-frame power-law diagnostic diagram; see §2.3 
and Figure 4. The tracks on this diagram show the ap- 
proximate regions where typical AGNs are likely to lie; 
since we investigate rest-frame energies up to 20 keV, CX- 
OHDFN J 123636.7+621 156 with a redshift of only z = 0.555 
cannot be included in these analyses. The AGN-classified 
SMGs are clearly well represented by this diagram and, on the 
basis of our model tracks, 12 (« 80%) are heavily obscured 
(TVh ^ 10 23 cm -2 ). Under the assumption that the different X- 
ray spectral properties of the SMGs are due to absorption, the 
lower energy X-ray emission should be more strongly atten- 
uated than the higher energy X-ray emission. This is demon- 
strated in Figure 5 where the implied attenuations from the lu- 
minosity ratios in the 2-10 and 5-20 keV bands are in agree- 
ment with those expected from the rest-frame 5-20 keV spec- 
tral slopes and the results from Figure 4. 

The constraints on the absorption properties of individual 
sources are comparatively poor. To provide tighter overall 
constraints we performed joint spectral fitting of all of the 
sources in the three different obscuration classes indicated 
in Figure 5 (N H < 10 23 cm -2 , N H = 1-5 xlO 23 cm -2 , and 
iV H > 5 x 10 23 cm -2 ). We jointly fitted the rest-frame 2- 
10 keV and 5-20 keV spectral slopes of the sources, leav- 
ing the power-law normalisations as a free parameter for each 



FIG. 5. — Rest-frame 5-20 keV spectral slope versus rest-frame X-ray lu- 
minosity ratio ( ^ 2 ~ IQkcV ). The solid line indicates the track for our adopted 

AGN model with T= 1.8 (individual values of Nn are indicated). The 
dashed lines indicate the divisions of the three obscuration classes used in 
the subsequent analyses (N}i < 10 23 cm -2 , Nn = 1 — 5 x 10 23 cm~ 2 , and 
Nu > 5 x 10 23 cm -2 ). The implied absorbing column densities of the AGNs 
in the SMGs are in good agreement with Figure 4. 



source. The results of this joint spectral fitting are shown 
in Figure 6. The average absorbing column densities for the 
sources in each obscuration class are in excellent agreement 
with those implied from the individual X-ray spectral analy- 
ses. Furthermore, the average intrinsic X-ray spectral slopes 
estimated from the model tracks (T w 1.8) are in good agree- 
ment with those found for typical AGNs in the local Universe 
(e.g., Nandra & Pounds 1994; George etal. 2000). 

3.4. Rest-frame 2-20 keV spectra 

The previous analyses have provided the broad-band X-ray 
spectral properties of the AGN-classified SMGs but have not 
been sensitive to discrete emission features (e.g., Fe Koc emis- 
sion). Unambiguous Fe Ka emission has not been identified 
in any of the individual X-ray spectra and, due to poor pho- 
ton statistics, it is difficult to constrain the X-ray continuum 
either side of the emission line, which significantly affects the 
accuracy of Fe Ka constraints. It is not possible to improve 
greatly the Fe Ka constraints via joint spectral fitting since the 
critical parameters are the emission-line and X-ray continuum 
normalisations, which vary from source to source. However, 
we can improve the overall signal-to-noise ratio of the spectra 
by stacking the data for the sources in each obscuration class 
and searching for the direct signature of Fe Ka emission in 
these composite X-ray spectra. 

When combining different X-ray spectra it is important to 
take into account the effective area and effective exposure of 
each source. We used XSPEC to output the unbinned X-ray 
spectrum of each source, folded by its effective area and ef- 
fective exposure. We then redshift corrected each X-ray spec- 
trum to create a rest-frame spectrum for each source. These 
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FIG. 6. — Rest-frame 5-20 keV spectral slope versus rest-frame 2-10 keV 
spectral slope, showing the results of the joint X-ray spectral fitting of the 
sources in the three obscuration classes from Figure 5. The model tracks are 
the same as those shown in Figure 4. The average properties of the sources 
in the three obscuration classes validate the X-ray spectral analyses for indi- 
vidual sources. 



rest-frame spectra were combined and then binned to increase 
the signal-to-noise ratio. Since the spectral response of Chan- 
dra ACIS-I varies with energy, in principal it is necessary 
to take into account the spectral response for each source 
when their spectra are redshifted. However, in the case of 
our sources the limiting factor is photon statistics rather than 
spectral resolution (e.g., see §3.3 of Alexander et al. 2003b). 

Composite rest-frame 2-20 keV spectra for the three differ- 
ent obscuration classes are shown in Figure 7. The signal-to- 
noise ratio of these X-ray spectra is significantly higher than 
the individual X-ray spectra and discrete features and the ac- 
curate shape of the X-ray continua can be seen. The obscured 
AGNs (Nr > 10 23 cm ) show a clear deficit of emission at 
< 5 keV due to absorption. As expected, there is a good 
correspondence between the amount of absorption and obscu- 
ration class. We plotted the model spectra determined from 
the joint X-ray spectral fitting (see Figure 6 and §3.3) and 
normalised them to the composite X-ray spectra by eye; see 
Figure 7. Reassuringly, the model tracks provide a good de- 
scription of the composite X-ray spectra, validating the results 
from the X-ray spectral analyses (Figures 3-6) and showing 
that the variations in the properties of the individual X-ray 
spectra are due to absorption. 

An emission feature at the rest-frame energy of Fe Koc 
is seen in the composite spectrum of the most heavily ob 



scured AGNs (Nn > 5 x 10 23 cm z ); see Figure 7. The rest- 
frame equivalent width of this feature (« 1 keV) is consis- 
tent with that expected for Compton-thick or near Compton- 
thick AGNs (e.g., Bassani et al. 1999) and suggests that some 
of the X-ray emission is reflected and/or scattered (e.g., Matt, 
Brandt, & Fabian 1996). The rest-frame energy of this fea- 
ture is closer to 6.7 keV than 6.4 keV, possibly indicating that 
it is produced in a warm scattering medium rather than from 
cold reflection, or that contributions from both are present. 



However, given the poor signal-to-noise ratio and effective 
spectral resolution of the data, we cannot distinguish between 
these possibilities (e.g., compare the different signal-to-noise 
ratio spectra of NGC 1068 and NGC 2992 in Figure 4 of 
Turner etal. 1997a). If all of the X-ray emission was re- 
flected/scattered then we would not be able to get an accu- 
rate estimate of the intrinsic X-ray luminosity of the AGNs 
in this obscuration class. However, at least for the properties 
of our model, a significant fraction of the direct X-ray emis- 
sion is seen at 10 keV, providing a reasonable constraint on 
the underlying AGN luminosity. There is ambiguous evidence 
for Fe Ka in the composite spectra of the Nn < 5 x 10 23 cm~ 2 
sources, potentially indicating the presence of Compton-thick 
or near Compton-thick AGN emission that is additionally 
obscured by starburst regions (e.g., Fabian etal. 1998). Im- 
proved photon statistics are required to test this scenario; see 
§4.4. 

In Figure 7 there appears to be an excess of < 4 keV emis- 
sion with respect to the model in the composite spectrum of 
the most heavily obscured sources (JV H > 5 x 10 23 cm" 2 ), 

:0.5- 



which may be due to star formation. The extrapolated 
8.0 keV luminosity (« 10 42 erg s -1 ) is consistent with the 
X-ray properties of the SMGs individually classified as star- 
burst galaxies; see Table 1 and §3.2. We further investigate 
the X-ray emission from star formation in §4.2. 

3.5. Absorption Correction 

Since accurate absorption corrections are challenging to de- 
termine even for well-studied obscured AGNs in the local 
Universe (e.g., Turner etal. 1997b; Bassani etal. 1999; Matt 
etal. 2000), the absorption corrections for our X-ray faint 
SMGs will be somewhat uncertain. Thankfully, the compos- 
ite X-ray spectra have shown that our AGN model gives a 
good characterisation of the X-ray properties of the AGN- 
classified SMGs and therefore provides a good base for de- 
termining plausible absorption corrections. In Table 2 we 
show our estimated unabsorbed X-ray luminosities for each 
of the AGN-classified SMGs using our AGN model and the 
column-density constraints from Figures 4-5. Although the 
absorption corrections for the most heavily obscured objects 
(JVh > 5 x 10 23 cm -2 ) are the most uncertain, they are consis- 
tent with those previously estimated for similar sources with 
comparable column densities and reflection/scattering com- 
ponents (e.g., Iwasawa etal. 2001, 2005b; Fabian etal. 2003; 
Wilman et al. 2003). Perhaps reassuringly, the mean unab- 
sorbed luminosities of the Nu > 5 x 10 cm -2 sources are 
consistent with those found for the Nh = 1-5 x 10 23 cm -2 
sources (L\ ~6 x 10 43 erg s versus L\ ~ 5 x 10 43 erg s _1 ), 
as expected if the differing amounts of absorption are primar- 
ily due to the orientation of the obscuration with respect to 
the central source (i.e., the unified AGN model; Antonucci 
1993). To provide a comparison to further studies in the lit- 
erature (e.g., Bautz etal. 2000; Mainieri etal. 2002; Szokoly 
et al. 2004; Severgnini et al. 2005), we also calculated unab- 
sorbed X-ray luminosities under the assumption of a simple 
absorbed power-law model (i.e., without the reflection and 
scattering components); see Table 2. However, we note that 
this model is unable to explain the presence of Fe Ka emis- 
sion in the most heavily obscured AGNs. 

Although our X-ray observations probe rest-frame ener- 
gies ^ 20 keV, it is possible that a large fraction of 
the X-ray emission is still hidden (e.g., much of the ob- 
served emission could be scattered rather than direct). More 
sensitive X-ray observations, and observations at > 10- 
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FIG. 7. — Composite rest-frame 2-20 keV spectra for each obscuration class, as indicated; see §3.3. The spectra are normalised to the average flux density 
of each obscuration class (with the exception of the the /% < 10 23 cm~ 2 sources, which have been scaled by a factor of 2 for presentation purposes) and are 
binned at 20 counts per bin. The total number of counts for each obscuration class are ks 690 (JV H < 10 23 cm" 2 ), ss 990 (Nr = 1-5 x 10 23 cm~ 2 ), and ss 580 
(Nh > 5 x 10 23 cm -2 ). The different line styles show the X-ray spectra for our adopted model for different amounts of X-ray absorption; the absorbing column 
density is taken from the joint X-ray spectral fitting for each obscuration class (see Figure 6). The Fe Ka line is at 6.4 keV in the model spectra; the apparent 
emission feature bluewards of the Fe Ka line is due to the supposition of the different continuum components and the Fe absorption edge at « 7.1 keV. The 
apparent excess of < 4 keV emission, with respect to the model, in the composite spectrum of the most heavily obscured sources (Af H > 5 x 10 23 cm" 2 ) may be 
due to star formation (see §3.4). 



20 keV (observed-frame energies), would be able to test 
this scenario (e.g., see Franceschini et al. 2000; Iwasawa 
etal. 2001 for the brighter but similarly luminous galaxy 
IRAS 09104+4109). Significant high-energy constraints 
could become available with NuSTAR as soon as 2009, be- 
fore further constraints with Constellation-X and XEUS in 
> 10 yrs time (see §4.4). 6 An alternative approach is to 
seek independent estimates of the AGN luminosity at other 
wavelengths (e.g., via the [OIII]A.5007 emission line or any 
strong near-IR to mid-IR high-excitation lines; Mulchaey 

6 On the basis of the expected ss 3 o" sensitivity limits for an ultra-deep 
» 1 Ms exposure (10-40 keV fluxes of ss 10~ 14 erg cm~ 2 s ; D. Stern, pri- 
vate communication), NuSTAR should detect emission down to Lx ~ 10 44 - 
10 45 ergs _I al z~ 1-3 ; see http://www.nustar.caltech.edu/ for further details. 



etal. 1994; Bassani etal. 1999; Sturm etal. 2002). Only one 
of our AGN-classified SMGs has an independent AGN lu- 
minosity constraint ([OIII] emission is detected from CXO- 
HDFN J123549.4+621536; T. Takata, in preparation). The 
predicted rest-frame 0.5-8.0 keV luminosity of this source 
(L x w 0.1-10 x 10 44 erg s _1 , determined using the [OIII]-X- 
ray relationship given in Bassani et al. 1999) is consistent with 
our estimate (Lx ~ 10 44 erg s _1 ); however, there are large 
uncertainties and this estimate is likely to increase when the 
effect of host-galaxy absorption towards the [OIII] emission- 
line region is taken into account. Further independent esti- 
mates of the AGN luminosity in SMGs would provide tighter 
constraints on the amount of nuclear absorption; see also §4.3. 



THE X-RAY SPECTRAL PROPERTIES OF SCUBA GALAXIES 



9 



4. DISCUSSION 

The analyses in §3 have provided a comprehensive char- 
acterisation of the X-ray properties of radio-selected spectro- 
scopically identified SMGs with fsso/na ~ 4 mJy. These prop- 
erties can be used to explore the role of AGN activity in these 
galaxies. In §4.1 we further investigate the AGN properties 
of our SMG sample, in §4.2 we estimate the contribution that 
AGN activity makes to the bolometric luminosity, in §4.3 we 
consider the constraints on black-hole growth, and in §4.4 we 
discuss the prospects for deeper X-ray observations and ex- 
plore the scientific potential of the next generation of X-ray 
observatories for the study of SMGs. 

4.1. AGN Activity in the SCUBA Galaxy Sample 

The X-ray properties of the AGN-classified SMGs in our 
sample are generally consistent with those of nearby lumi- 
nous AGNs (i.e., T « 1.8 ±0.5, N H « 10 20 -10 24 cm" 2 , and 
L x « 10 43 -10 44 5 erg s _1 ; e.g., Smith & Done 1996; Turner 
etal. 1997b; Maiolino etal. 1998; Risaliti etal. 1999); see 
§3 and Table 2. The majority (« 80%) of the AGNs are 
heavily obscured (Nh ~ 10 23 cm -2 ), four of which have 
L x ~ 10 44 erg s -1 and could be considered obscured (Type 2) 
quasars; see Bautz et al. (2000) and Mainieri et al. (2005a) for 
other examples of obscured quasars that are bright at submm 
wavelengths. With the exception of completely Compton- 
thick AGNs (i.e., A'h ~ 10 25 cm -2 ), the overall column den- 
sity distribution would be roughly similar to that found for 
nearby AGNs (e.g., Risaliti etal. 1999; Salvati & Maiolino 
2000). The lack of obvious X-ray emission from AGN ac- 
tivity in the X-ray classified starbursts could be because the 
AGNs are of a low luminosity (Lo.5-8.0keV ~ 10 42 erg s _1 ) 
or are completely Compton thick (A^h « 10 25 cm~ 2 ). 7 In- 
terestingly, if the latter is assumed then the estimated column 
density distribution would be remarkably similar to that found 
for nearby AGNs. 

The obscured AGN fraction in our sample (« 80% have 
A^h £ 10 23 cm~ 2 ) is larger than that found from identifica- 
tion studies of the w 2-8 keV background population (e.g., 
Ueda et al. 2003 find similar numbers of obscured and unob- 
scured AGNs; see also Barger et al. 2002; Mainieri et al. 2002; 
Szokoly etal. 2004). Although this may suggest a prevalence 
for obscured activity in the SMG population, selection effects 
are also likely to be at least partially responsible. For exam- 
ple, due to the comparatively high redshifts of our sources we 
probe considerably higher rest-frame energies than those for 
the bulk of the spectroscopically identified sources in deep 
X-ray surveys [e.g., « 1.5-24 keV (at z ~ 2) as compared 
to w 0.9-14 keV (at z ~ 0.7) at observed energies of 0.5- 
8.0 keV], allowing higher column densities to be penetrated. 
Further z > 1 obscured AGNs are likely to have been de- 
tected in X-ray surveys but lack spectroscopic redshifts due 
to the faintness of their optical counterparts (e.g., Alexan- 
der etal. 2001; Barger etal. 2002; Fiore etal. 2003; Treister 
etal. 2004; Mainieri etal. 2005b). These optically faint X- 
ray sources are the dominant population of obscured AGNs 
at z > 1 and may account for as much as w 50% of the AGNs 
detected in deep X-ray surveys (Alexander et al. 2002). Com- 
plete spectroscopic identification of these sources is required 
to determine the obscured to unobscured AGN ratio at z > 1 

7 Of course some of the X-ray classified starbursts may not host AGN 
activity. However, we note that the similar z = 2.285 galaxy FSC 10214+4724 
hosts an AGN that is weak at X-ray energies and would have been classified 
as a starburst galaxy in our study (Alexander etal. 2005b). 



and reveal whether our X-ray detected SMGs are preferen- 
tially more obscured than the general z > 1 X-ray source pop- 
ulation. 

4.2. What Powers SCUBA Galaxies? 

AGN activity clearly plays an important role in the spec- 
troscopically identified SMG population but does it domi- 
nate the bolometric output? The answer to this key ques- 
tion can have important implications for the formation and 
evolution of massive galaxies and the growth of massive 
black holes (see §1). As found for most dusty luminous 
galaxies, the luminosity of SMGs is likely to peak at far-IR 
wavelengths and have components of both AGN and star- 
formation activity (e.g., Sanders & Mirabel 1996; Genzel & 
Cesarsky 2000). Many studies have suggested that star for- 
mation accounts for a large fraction of the bolometric out- 
put of SMGs (e.g., Frayer et al. 1998, 2004; Ivison et al. 2002, 
2004; Alexander etal. 2003b; Almaini etal. 2003; Chapman 
etal. 2004a; Egami et al. 2004; Swinbank et al. 2004) but none 
has resolved the rest-frame far-IR emission. The most direct 
study to date is that of Chapman et al. (2004a) who showed 
that w 70% of radio-identified SMGs have radio emission ex- 
tended on w 10 kpc scales, implying that the rest-frame far-IR 
emission is also extended over these scales. However, these 
results do not rule out the presence of energetically significant 
AGN components on smaller scales. The important contribu- 
tion our study can make to this debate is in providing direct 
estimates of the bolometric contribution from AGN activity. 

In Figure 8 we show the rest-frame far-IR luminosity ver- 
sus unabsorbed 0.5-8.0 keV luminosity for the SMGs and 
compare them to well-studied starburst galaxies and AGNs 
drawn from the literature. The shaded region indicates the 
typical range of luminosity ratios for the well-studied quasars 
of Elvis et al. (1994) and provides an indication of the location 
of AGN-dominated sources on this figure. None of the far-IR 
luminous SMGs (Lfir ^ 10 12 L ) lies within the dark-shaded 
region, indicating that they are comparatively weak at X-ray 
energies. The difference between the median X-ray-to-far-IR 
luminosity ratios of the AGN-classified SMGs 0.004) 

and the quasars (j^~ 0.05) suggest that the AGN activity in 
the SMGs contributes, on average, only « 8% of the far-IR 
emission. This comparison indicates that star formation typi- 
cally dominates the bolometric output of these SMGs. How- 
ever, if the SMGs have an AGN dust-covering factor « 12 
times larger than found in the quasars then they could poten- 
tially be AGN dominated (i.e., a larger amount of dust could 
be heated for a given central source luminosity). Although 
the large fraction of obscured AGNs in our sample suggests 
that the dust-covering factor is large in these SMGs, the dust- 
covering factor of quasars is uncertain and may depend on the 
evolutionary state or nature of individual objects (e.g., Haas 
et al. 2003; Page et al. 2004). 

We can take a different approach and compare our SMGs to 
the literature galaxies, many of which host an obscured AGN 
and may be physically similar to the SMGs; see Figure 8. The 
complication of this comparison is that the relative contribu- 
tion from AGN and star-formation activity in these galaxies is 
often poorly constrained. Even so, we can gain some insight 
into the relative dominance of AGN and star-formation activ- 
ity by "calibrating" Figure 8 using the ISO mid-IR spectral 
diagnostics of Rigopoulou et al. (1999) and Tran et al. (2001); 
see also Genzel etal. (1999) and Lutz, Veilleux, & Genzel 
(1999). On the basis of these studies, we have found that the 
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FIG. 8. — Rest-frame far-IR versus unabsorbed 0.5-8.0 keV luminosity for our SMGs (circles) and galaxies selected from the literature (squares; individual 
galaxies are labelled to allow objects to be compared to the SMGs). Sources are further highlighted as follows: AGN-classified SMGs (filled circles), starburst- 
classified SMGs (open circles), literature galaxies (filled squares); literature galaxies classified as AGN-dominated or star-formation dominated by Rigopoulou 
et al. (1999) and Tran et al. (2001) are additionally indicated by an "A" or "S", respectively. The diagonal lines show ratios of constant X-ray-to-far-IR luminosity: 
the dotted line shows the mean luminosity ratio of the starburst-classified SMGs, the solid line shows the median luminosity ratio of the AGN-classified SMGs, 
and the dashed line shows the median luminosity ratio for the quasars studied by Elvis etal. (1994). The shaded region indicates the standard-deviation in 
the luminosity ratio of the quasars studied by Elvis etal. (1994). The X-ray and infrared data for the literature galaxies were taken from Sanders & Mirabel 
(1996), Hughes etal. (1997), Bassani etal. (1999), Bautz etal. (2000), Matt etal. (2000), Risaliti etal. (2000), Iwasawa (2001), Iwasawa, Fabian, & Ettori (2001), 
Gallagher etal. (2002), Lira etal. (2002), Verma etal. (2002), Braito etal. (2003, 2004), Franceschini etal. (2003), Sanders etal. (2003), Alexander etal. (2005b), 
Brandt & Hasinger (2005), Iwasawa etal. (2005b), and the NASA Extragalactic Database (NED). The literature data have been converted to the cosmology used 
here and, when appropriate, the X-ray data have been converted to the 0.5-8.0 keV band assuming V = 1.8 for the AGNs and T = 2.0 for the starburst galaxies. 
Although the unabsorbed X-ray luminosities are presented whenever possible, in some sources the intrinsic power of the AGN is unknown (e.g., the possible 
Compton-thick AGN in FSC 10214+4724; Alexander etal. 2005b) and the observed X-ray luminosity is shown. In the absence of far-IR data, far-IR luminosities 
were estimated from the 1 .4 GHz luminosity density under the assumption of the radio-to-far- IR relationship. 



threshold between the AGN-dominated and star-formation 
dominated literature galaxies is j^L ~ 0.004; of the 35 galax- 
ies in Rigopoulou et al. (1999) and Tran et al. (2001), all of the 
17 AGN-dominated galaxies have ^ > 0.004 while 13 of 

the 18 star-formation dominated galaxies have j^- < 0.004 
(see Figure 8). This indicates that AGN-dominated galaxies 
can lie outside of the shaded region in Figure 8 and suggests 
that the Elvis et al. (1994) bolometric correction is unlikely to 
be applicable for all sources. The implication of these results 



is that some of the SMGs could be predominantly powered by 
AGN activity even though the majority are likely to be star- 
formation dominated. However, since different techniques of- 
ten yield contrasting results (e.g., Farrah et al. 2003), and the 
exact contributions from AGN and star-formation activity are 
uncertain even in extremely well-studied local galaxies (e.g., 
Spoon et al. 2000; Lutz et al. 2003; Iwasawa etal. 2005a), this 
conclusion can only be considered qualitative. Deep mid-IR 
spectroscopy with the Spitzer telescope should provide im- 
portant additional insight into the relative contributions from 
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AGN and star-formation activity in these SMGs (e.g., Lutz 
et al. 2005 for the first constraints). However, significantly im- 
proved constraints are unlikely to be achieved before the ad- 
vent of the Atacama Large Millimeter Array (ALMA), which 
will offer spatial resolution constraints at submm wavelengths 
on w 100 pc scales for galaxies at z ~ 2. 8 

In Figure 8 we also show the X-ray and far-IR luminosi- 
ties of the starburst-classified SMGs. The mean X-ray-to- 
far-IR luminosity ratio of these sources (jj^ ~ 10~ 4 , when 
the X-ray stacking result of the three X-ray undetected SMGs 
is included; see §2.2) is consistent with that expected from 
star formation if the X-ray emission is dominated by HMXBs 
(Persic etal. 2004; see also Figure 2b and §3.2). Under this 
assumption the X-ray derived star-formation rates (m 1300- 
2700 M yr~') are in good agreement with those determined 
using other techniques (e.g., Swinbank etal. 2004; Chapman 
et al. 2005). The mean X-ray-to-far-IR luminosity ratio is w 4 
times lower than that found for typical starburst galaxies in the 
local Universe. However, as argued in Persic et al. (2004; see 
also Franceschini et al. 2003), the X-ray emission from typical 
starburst galaxies will be dominated by long-lived LMXBs, 
which will be comparatively weak in intense star-forming 
galaxies. Of course, the weak X-ray emission from these 
sources could also imply that luminous completely Compton- 
thick AGNs are present (e.g., Iwasawa et al. 2005a). 

4.3. Growth of Black Holes in Massive Star-Forming 
Galaxies 

In Alexander etal. (2005a) we used the results presented 
here to constrain the growth of black holes in massive star- 
forming galaxies. We showed that the AGN fraction in 
the bright SMG population is > 38+}q%, when corrected 
for bright SMGs without spectroscopic redshifts, and argued 
that this indicates that their black holes are almost continu- 
ously growing throughout vigorous star-formation episodes. 
The most likely catalyst for this activity appeared to be 
galaxy major mergers (e.g., Chapman etal. 2003c; Conselice 
et al. 2003). Shortly after the completion of our study sophis- 
ticated hydrodynamical simulations of galaxy major mergers, 
taking into account the growth of both the black hole and stel- 
lar components, were published (e.g., Di Matteo etal. 2005; 
Springel et al. 2005). The simulation results are in good agree- 
ment with the conclusions of Alexander et al. (2005a), with 
the peak epoch of star formation corresponding to a heav- 
ily obscured rapid black-hole growth phase, which is ul- 
timately proceeded by an unobscured quasar phase (Hop- 
kins etal. 2005). However, the black holes of the most 
massive galaxies in these simulations (which are probably 
most similar to the SMGs) are up-to an order of magni- 
tude more massive than those estimated from our X-ray lu- 
minosities under the assumption of Eddington-limited accre- 
tion (M BH ( Edd ) 10 8 M ). Our results can be better rec- 
onciled with the simulations if either the black holes in the 
AGN-classified SMGs are accreting at sub-Eddington rates, 
a large fraction of the AGN activity is hidden at rest-frame 
energies of ss 20 keV (see §3.5), or the X-ray-to-bolometric 
luminosity ratio is significantly lower than that typically as- 
sumed for AGNs (i.e., Elvis et al. 1994; see Table 1 of Alexan- 
der etal. 2005a for further details). However, there is further 
observational evidence that the masses of the black holes in 
SMGs may be comparatively modest and in good agreement 



See http://www.alma.nrao.edu/ for more details on ALMA. 



with our estimate. Rest-frame optical spectra of SMGs show 
that when broad emission lines are detected they are often 
comparatively narrow (typical full width at half-maximum 
velocities of w 1000-3000 km -1 ; Vernet & Cimatti 2001; 
Smail et al. 2003; Swinbank et al. 2004, 2005; T. Takata et al. , 
in preparation). Under the assumption that the dynamics of 
the broad-line regions in these sources are dominated by the 
gravity of the central black hole, these comparatively nar- 
row emission line widths suggest typical black-hole masses of 
10 8 M for mass accretion rates of 1 M yr~' (Alexan- 
der et al. 2005a; e.g., McLure & Dunlop 2004). Although it is 
premature to draw conclusions before more accurate black- 
hole mass constraints are secured, we note that other models 
have suggested that the masses of the black holes in SMGs are 
comparatively modest (e.g., Archibald et al. 2002; Kawakatu 
etal. 2003; Granato etal. 2004a,b). In particular, the physi- 
cally motivated models of Granato et al. bear a striking simi- 
larity to our results (see Figures 2 & 3 in Granato et al. 2004b 
for the M vir = 10 12 4 M model). 

4.4. Prospects for Deeper X-ray Observations 

The quality of the individual X-ray spectra are restricted by 
poor photon statistics. However, the composite X-ray spec- 
tra in §3.4 demonstrate what could be achieved for individ- 
ual sources with extremely long Chandra exposures (up to 
k 12 Ms). With individual X-ray spectra of this quality we 
would be able to significantly improve the constraints for in- 
dividual sources (e.g., F, Nh, the strength of Fe Koc emis- 
sion, intrinsic X-ray luminosity, measure soft X-ray excesses 
due to star-formation activity). The increased photon statis- 
tics would also allow for basic X-ray spectral analyses (up to 
k 200 counts) of the starburst-classified SMGs. 

The major scientific thrust of the next generation of X-ray 
observatories {Constellation-X and XEUS) is high-resolution 
spectroscopy of faint X-ray sources. 9 In order to investigate 
the scientific potential of these observatories for the study of 
SMGs, we have simulated a 300 ks XEUS exposure of one of 
the most heavily obscured sources (SMMJ 123622.6+621629 
at z = 2.466) using XSPEC; see Figure 9. An Fe Koc emission 
line is easily identifiable at rest-frame 6.4 keV, indicating that 
this is a Compton-thick or near Compton-thick AGN; similar 
results can be obtained with Constellation-X for sources « 10 
times brighter. The quality (signal-to-noise ratio and spectral 
resolution) of the simulated X-ray spectrum exceeds that of 
our composite X-ray spectra (compare with Figure 7). Both 
Constellation-X and XEUS will also have the capability to de- 
tect sources at > 10-20 keV (observed-frame energies) and 
place crucial AGN luminosity constraints at very high X-ray 
energies (see Footnote 6 for NuSTAR constraints); see §3.5. 

The increased photon statistics and spectral resolution of 
the next-generation X-ray observatories should allow for the 
detection of weaker AGN features. For example, some lo- 
cal AGNs and powerful quasars have shown evidence for 
large-scale outflowing material, probably due to accretion- 
disk winds (e.g., Chartas etal. 2002; Kaspi etal. 2002; Ogle 
etal. 2003). These outflows provide an efficient method for 
distributing high-metallicity gas into the intergalactic medium 
and may be responsible for producing the M-o relationship 
seen in local galaxies (e.g., Silk & Rees 1998; Fabian 1999; 
King 2003; Di Matteo et al. 2005). If SMGs are the precursors 
to luminous quasar activity, and the ancestors of nearby 

9 See http://constellation.gsfc.nasa.gov/science/index.html and 
http://www.rssd.esa.int/XHUS/ for more details. 
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Observed Energy (keV) 

FIG. 9. — Simulated 300 ks final-configuration XEUS spectrum of the heavily obscured source SMMJ 123622.6+621629 (top panel). The spectrum was 
simulated using our adopted AGN model (see §2.3) with T = 1.8 and Nu = 1.5 x 10 24 cirT 2 , and additionally includes an accretion-disk wind outflow [taken 
from NGC 1068 (Ogle etal. 2003) and scaled to the Fe Ka flux of our source]; see §4.4. The solid line indicates the input model with the emission-line 
components removed to highlight the emission-line features. The bottom panel shows the ratio of the simulated data to the model. Rest-frame 6.4 keV Fe Ka 
is easily identified (at sa 1.8 keV) and indicates that this source is a Compton-thick or near Compton-thick AGN. The < 1 keV emission features are due to 
the outflowing accretion-disk wind. Six of the strongest emission lines are labelled in the bottom panel (see Ogle et al. 2003 for further information and line 
identifications). 



massive galaxies, then they are likely to produce signifi- 
cant emission-line outflows (see Smail etal. 2003; Chapman 
et al. 2004b; S winbank et al. 2005 for evidence at optical/near- 
IR wavelengths). We included the emission-line outflow pa- 
rameters found for NGC 1068 (Ogle etal. 2003) to our sim- 
ulated 300 ks XEUS spectrum, scaled to the Fe Ka flux of 
our source; see Figure 9. The emission-line signatures of the 
outflow are identifiable in the simulated spectrum at £ 1 keV, 
allowing basic constraints to be placed on the properties of the 
outflowing material (e.g., the mass outflow rate). These con- 
straints could prove to be crucial in our understanding of the 
connection between AGN and star-formation activity at high 
redshift, and help to explain the properties of nearby massive 
galaxies. 

5. CONCLUSIONS 

Using the powerful combination of ultra-deep X-ray obser- 
vations and deep optical spectroscopic data we have placed 
constraints on the X-ray properties of 20 radio-selected spec- 
troscopically identified SMGs with fsso/jm ~ 4 mJy. Our key 
results are the following: 

1 . Seventeen of the 20 SMGs in our sample are detected 
at X-ray energies. From a classification of their X-ray 
properties, we found the X-ray emission to be domi- 
nated by AGN activity in 15 sources (~ 75%); the X- 
ray emission from the other 5 sources is likely to be 
dominated by star-formation activity. See §3.2. 

2. Using a variety of X-ray spectral analyses we 
found that the properties of the AGNs are gener- 
ally consistent with those of nearby luminous AGNs 



(i.e., T « 1.8 ±0.5, N H « 10 20 -10 24 cm" 2 , and 
L x « 10 43 -10 44 5 erg s" 1 ); see §3 and Table 2. 
The majority (~ 80%) of the AGNs are heavily 
obscured (Nh £ 10 23 cm -2 ), four of which have 
Lx ~ 10 44 erg s _1 and could be considered obscured 
quasars. The estimated column-density distribution is 
roughly similar to that found for nearby AGNs. See 
§3.3 and §4.1. 

3. We constructed rest-frame 2-20 keV composite spec- 
tra for three different obscuration classes (A^h < 
10 23 cm" 2 , N H = 1-5 xlO 23 cm" 2 , and N H > 5 x 
10 23 cm -2 ). These composite spectra revealed fea- 
tures not observed in the individual X-ray spectra and 
an w 1 keV equivalent width Fe Ka emission line is 
seen in the composite X-ray spectrum of the most heav- 
ily obscured AGNs, suggesting Compton-thick or near 
Compton-thick absorption. The good agreement be- 
tween our model tracks and the composite X-ray spec- 
tra validate the X-ray spectral analyses of individual 
sources. See §3.4. 

4. Taking into account the effects of absorption, we found 
that the average X-ray to far-IR luminosity ratio of 
the AGN-classified sources (j^= 0.004) is approxi- 
mately one order of magnitude below that found for 
typical quasars. This result suggests that intense star- 
formation activity (of order « 1000 M©yr _1 ) domi- 
nates the bolometric output of spectroscopically iden- 
tified SMGs. Possible biases in the selection of this 
sample (see §2.1) would suggest that, if anything, they 
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should have a higher AGN contribution than spectro- 
scopically unidentified SMGs. We also investigated the 
possibility that the X-ray to far-IR luminosity ratio for 
the AGNs in SMGs is intrinsically less than that found 
for typical quasars and postulated that some SMGs may 
be AGN dominated. See §4.2. 

5. We found that the X-ray emission from the starburst- 
classified SMGs is consistent with that expected from 
HMXBs. The X-ray derived star-formation rates 
(~ 1300-2700 M©yr _1 ) are in good agreement with 
those determined using other techniques. See §4.2. 

6. We find good agreement between our overall picture for 
the growth of black holes in massive galaxies and the 
results from recent hydrodynamic simulations. We pro- 
vide support that the black holes in our AGN-classified 
SMGs are < 10 8 M . See §4.3. 

7. We demonstrated that the next generation of X-ray ob- 
servatories (Constellation-X and XEUS) will place sig- 
nificantly improved X-ray spectral constraints on in- 
dividual sources. In addition to constraining the ab- 
sorption properties of the AGN components we found 
that these X-ray observatories have the potential to de- 



tect the presence and properties of outflowing material. 
These observations could prove to be crucial in our un- 
derstanding of the evolution of massive galaxies. See 
§4.4. 
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Table 1 . Basic properties of the SCUBA galaxies 
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" SCUBA galaxy properties, taken from Chapman et al. (2005). Submm flux density is in units of mjy. 

b Rest-frame 1.4 GHz luminosity density in logarithmic units of W Hz -1 . Calculated following equation 2 in Alexander et al. (2003b) for a = 0.8 (the average spectral index for star-forming galaxies; e.g., Yun, 
Reddy, & Condon 2001) using the radio flux density from Chapman et al. (2005). 

c Rest-frame far-IR (k =40-120 /jm) luminosity in logarithmic units of solar luminosities. Calculated from the rest-frame 1.4 GHz luminosity under the assumption of the radio to far-IR correlation with q = 2.35 
(e.g., Helou, Soifer, & Rowan-Robinson 1985). 

d X-ray source identification number and X-ray source co-ordinates; taken from Alexander et al. (2003a). "SUPP" refers to a source detected in the complete supplementary Chandra catalog; see Alexander et al. 
(2003b), §3.2 and §3.4.2 of Alexander et al. (2003a). 

c Ratio of the count rates in the 2.0-8.0 keV and 0.5-2.0 keV bands, and the effective photon index in the 0.5-8.0 keV band (calculated from the band ratio); taken from Alexander et al. (2003a). T = 1.4 has been 
adopted for sources with poorly defined band ratios. 

f Full-band flux uncorrected for Galactic or intrinsic absorption, in units of 10~' 5 erg cm -2 s _1 . Taken from Alexander et al. (2()03a,b). Upper limits are calculated following §3.4.1 of Alexander et al. (2003a). 
g Rest-frame 0.5-8.0 keV luminosities in logarithmic units of erg s -1 , calculated following equation 1 in Alexander et al. (2003b) for T = 1.8.. 

h Source classification based on X-ray properties; see §3.2. "AGN" indicates the X-ray emission is dominated by AGN activity, "(V)" indicates evidence for X-ray variability, and "S/Burst?" indicates that the X-ray 
emission is probably dominated by luminous star-formation activity. <-n 



Table 2. Spectral fits for the AGN-classified SCUBA galaxies 
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a X-ray source position, taken from Table 1 . 

b Corresponding rest-frame energy (keV) for the observed 0.5-10 keV band. 
c Total number of 0.5-10 keV counts used in the X-ray spectral analyses. 

d Best-fit photon index for a simple power-law model (with Galactic absorption) in the indicated rest-frame energy bands; the uncertainties refer to the 90% confidence level (for one interesting parameter). 

e Estimated X-ray absorbing column density in logarithmic units of cm -2 . See §3 for justification and column density analyses. 

f Rest-frame X-ray luminosity, determined directly from XSPEC using the best-fit photon indices. This has not been corrected for absorption. 

g Rest-frame 0.5-8.0 keV luminosity in logarithmic units of erg s _1 , determined from the rest-frame 5-20 keV luminosity and corrected for the estimated X-ray absorbing column density using our adopted AGN 
model; see §3.5. 

h Rest-frame 0.5-8.0 keV luminosity in logarithmic units of erg s _1 , determined from the rest-frame 5-20 keV luminosity and corrected for the estimated X-ray absorbing column density using a simple absorbed 
power-law model; see §3.5. 



